A Search for Electron Antineutrino Appearance at the Am 2 ~ 1 eV 2 Scale 
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(Dated: September 28, 2009) 

The MiniBooNE Collaboration reports initial results from a search for P M — > i? e oscillations. A 
signal-blind analysis was performed using a data sample corresponding to 3.39 x 10 20 protons on 
target. The data are consistent with background prediction across the full range of neutrino energy 
reconstructed assuming quasielastic scattering, 200 < E® E < 3000 MeV: 144 electron-like events 
have been observed in this energy range, compared to an expectation of 139.2 ± 17.6 events. No 
significant excess of events has been observed, both at low energy, 200-475 MeV, and at high energy, 
475-1250 MeV. The data are inconclusive with respect to antineutrino oscillations suggested by data 
from the Liquid Scintillator Neutrino Detector at Los Alamos National Laboratory. 



Motivated by the LSND observation of an excess of D e 
events in a beam [l[ , the MiniBooNE collaboration has 
previously performed a search for z/„ — > v e oscillations, 
the results of which showed no evidence of an excess of 
v e events for neutrino energies above 475 MeV @, Q ■ As- 
suming no CPT or CP violation, the results exclude the 
LSND excess interpreted as two-neutrino oscillations at 
Am 2 - 0.1-100 eV 2 at 98% C.L.. Similarly, the KAR- 
MEN experiment Q has performed a direct search for v e 
appearance, and has placed a limit independent of any 
CPT or CP violation assumption. However, a joint anal- 
ysis of KARMEN and LSND results shows high compat- 
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ibility 0, Q . A corresponding — > v e oscillation search 
has been performed at MiniBooNE and is presented in 
this Letter. This search serves as another direct test of 
LSND and provides complementary information to that 
of KARMEN, having sensitivity to the lower Am 2 oscil- 
lations allowed by the joint KARMEN-LSND analysis 0]. 
It should be noted that, in a simple two-neutrino oscilla- 
tion model where CPT conservation is imposed, oscilla- 
tion probabilities (mixing amplitudes and mass-squared 
differences) for neutrinos and antincutrinos cannot be dif- 
ferent. Therefore, the oscillation search presented here is 
a direct search for existence of non-standard oscillations 
where CPT is violated, or effectively violated. 

Despite having observed no evidence for oscillations 
above 475 MeV, the MiniBooNE — > v e search ob- 
served a 3.0cr excess of electron-like events at low energy, 
between 200-475 MeV Q • Although the excess is incom- 
patible with LSND-type oscillations, several hypotheses 
[l3| , including sterile neutrino oscillations 
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with CP violation, anomaly-mediated neutrino-photon 
coupling, and others, have been proposed that provide 
a possible explanation for the excess itself, and, in some 
cases, offer the possibility of reconciling the MiniBooNE 
v e excess with the LSND v e excess. These phenomenolog- 
ical interpretations have provided additional motivation 
for an antineutrino appearance search at MiniBooNE. 

The analysis presented in this Letter mirrors the blind 
search performed in neutrino mode 0j. It employs a 
two-neutrino oscillation model, where only present in 
the MiniBooNE beam are allowed to oscillate into v ei at 
Am 2 ~ 0.1-100 eV 2 . Given that no evidence of oscilla- 
tions was observed in high-purity, high-statistics searches 
in neutrino mode 0, EH ; the analysis further assumes no 
disappearance and no oscillations. In addition, no 
contribution from the observed neutrino mode low en- 
ergy excess has been accounted for in the antineutrino 
prediction. 

The antineutrino flux (IBj is produced by 8 GeV pro- 
tons incident on a beryllium target. Negatively charged 
mesons produced in p-Be interactions are focused in the 
forward direction with the use of a toroidal magnetic 
field, and subsequently decay primarily into v^. In an- 
tineutrino mode, a large neutrino contamination (y^ and 
v e ) of 15.9% is expected in the flux viewed by the detec- 
tor, compared to 5.9% in neutrino mode. The intrinsic v e 
and v e content is only 0.4% and 0.2%, respectively, com- 
ing primarily from 7r — > /i and K decays. The flux 
peaks at ~ 400 MeV and has a mean energy of ~ 600 
MeV. See [l5| for more details. 

A detailed description of the MiniBooNE detector is 
available in [l^]. The detector location was chosen to 
satisfy L[m]/E[MeV] ~ 1, similar to that of LSND, 
thus maximizing sensitivity to oscillations at Am 2 ~ 
1 eV 2 . The detector is filled with pure mineral oil 
(CH2). Neutrino interactions in the detector produce fi- 
nal state electrons or muons, which produce scintillation 
and Cherenkov light detected by photomultiplier tubes 
(PMTs) that line the interior of the detector. The sim- 
ulation of light incident on the PMTs takes into account 
decays and strong and electromagnetic re-interactions in 
the detector, and includes processes that were added in 
the final v e appearance analysis The M® E appear- 
ing in the nucleon axial vector form factor, and the Pauli 
blocking parameter, k, used to parametrize neutrino 
quasi-elastic scattering on carbon, were adjusted by fits 
to MiniBooNE data, as were the coherent pion cross sec- 
tions The m a E and K values of 1-23 ±0.08 GeV 



and 1.022±0.021, respectively, were used in this analysis. 
Two additional parameters, M^ E,H = 1.13 ± 0.10 GeV 
and Mp" = 1.10 ± 0.10 GeV, were introduced in the 
analysis to parameterize antineutrino quasi-elastic scat- 
tering on hydrogen and single pion production on hydro- 
gen. These processes have a non-negligible contribution 
in antineutrino running mode, where roughly 25% of the 
antineutrino quasi-elastic scatters are on hydrogen rather 



than carbon. 

The detector cannot differentiate (on an event-by- 
event basis) a from a interaction, or a v e from a v e 
interaction. Therefore, the reconstruction and selection 
requirements for £5 e -induced chargcd-currcnt quasi-elastic 
(CCQE) events, which is the characteristic signature of 



any possible signal from v„ 



oscillations, are identi- 



cal to those of the final neutrino mode analysis Q . 

To provide a constraint on v e candidate events, a P M 
CCQE sample is also formed by looking for events with a 
muon-like Cherenkov ring and a cluster of delayed PMT 
hits from the decay of the muon into an electron. The 
first cluster of PMT hits (muon subevent) is required to 
have more than 200 inner detector PMT hits, and no 
more than six outer (veto) PMT hits. A maximum of 
200 inner detector and six veto PMT hits are required 
for the second subevent (decay electron), and a mini- 
mum time cut of 1000 ns between the first and second 
subevents is required to ensure PMT stability for proper 
charge response. After reconstruction, the first subevent 
vertex and the track end-point under the muon hypothe- 
sis are required to occur within the fiducial volume. The 
neutrino energy reconstructed from the outgoing muon 
energy and angle, E® E , is required to satisfy E® E > 150 
MeV. A cut on the separation distance between the muon 
and decay electron vertices as a function of reconstructed 
energy of the muon is also applied to provide rejection 
against backgrounds, mostly from CC tt + interactions. 
For more details on the reconstruction method, see [lj)[ . 

The oscillation parameters are extracted from a com- 
bined fit to v e CCQE and CCQE event distributions, 
following This fit method takes advantage of strong 
flux and cross section correlations among the v e CCQE 
and CCQE event samples, since any possible — > v e 
signal, as well as some v e backgrounds, interact through 
the same process as CCQE events, and are related to 
CCQE events through the same n + or ir~ decay chain 
at production. These correlations enter through the off- 
diagonal elements of the covariance matrix used in the 
X 2 calculation, relating the contents of the bins of the 
v e CCQE and CCQE distributions. This procedure 



maximizes the sensitivity to 



oscillations when 



systematic uncertainties are included [201 ]. 

A sample of 14,107 data events passing P M CCQE se- 
lection requirements is used in the analysis. This sam- 
ple is compared to a MontcCarlo prediction which has 
been corrected to match the observed CCQE data 
through a normalization factor of 1.22 applied to events 
from 7T~ decays in the beam, and 0.93 applied to events 
from 7r + decays in the beam. These normalization fac- 
tors are extracted from a fit to the angular distributions 
of the outgoing /i + and fx~ in and CCQE interac- 
tions These two factors result in an overall 15% nor- 
malization correction which is covered by flux and cross 
section uncertainties. The same normalization correction 
is also applied to all possible signal events which share 
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TABLE I: The expected number of events for different E„ E 
ranges (in MeV) from all of the backgrounds in the i> e appear- 
ance analysis and for the LSND central expectation (0.26% 
oscillation probability) of — > v e oscillations, for 3.39 x 10 20 
POT. 
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The normal- 



ization correction is accounted for in the oscillation fit 
by a reduction in the quoted effective degrees of freedom 
by one unit. After correction, the sample contains 95% 
9ft and produced in pion decays and 2.4% 9^ and 
produced in kaon decays. The neutrino content of the 
sample is 22%. The majority of events (71%) are true 
CCQE interactions, with CCn^ 1 interactions being the 
dominant source of background (20%). This sample is 
included in the 9 e appearance fits as a function of 8 bins 
of reconstructed neutrino energy, E® E , ranging from 
to 1900 MeV. 

Table U shows the number of predicted 9 e CCQE back- 
ground events for different ranges of E® E . The back- 
ground estimates include both antineutrino and neu- 
trino events, the latter representing ~ 44% of the to- 
tal. The predicted backgrounds to the 9 e CCQE sam- 
ple are constrained by internal measurements at Mini- 
BooNE. These measurements use event samples from re- 
gions in reconstructed kinematic variables where any pos- 



sible signal from 



is negligible, in order to pre- 



serve blindness. The constrained backgrounds include 
NC 7T° events, A — > Nj radiative events, and events from 
interactions outside the detector. The NC w° background 
events arc adjusted in bins of 7r° momentum according 
to a direct tt° rate measurement in antineutrino mode, 
following (l7| . which uses events reconstructed near the 
7T° mass peak. The size of the applied correction to the 
total NC 7T° rate is less than 10%. The A -> V7 rate 
is indirectly constrained, being related to the measured 
7T° rate through a branching fraction and final state in- 
teraction correction. The rate of backgrounds from ex- 
ternal interactions is constrained through a direct mea- 
surement at MiniBooNE, using a separate event sample 
where the rate of external interaction events is enhanced. 
Other backgrounds from mis-identified or 9^ receive 
the 9 ft CCQE normalization correction according to their 
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FIG. 1: Top: The E® E distribution for u e CCQE data (points 
with statistical errors) and background (histogram with un- 
constrained systematic errors). Bottom: The event excess as 
a function of E® E . Also shown are the expectations from 
the best oscillation fit and from neutrino oscillation parame- 
ters in the LSND allowed region. The error bars include both 
statistical and systematic errors. 

parentage at production (ir + or ir~). Intrinsic v e and 9 e 
events from the 7r — > /i decay chain also receive this nor- 
malization. 

Systematic uncertainties are determined by consider- 
ing the effects on the 9^ and 9 e CCQE rate prediction 
of variations of fundamental parameters within their as- 
sociated uncertainty. These include uncertainties on the 
flux estimate, including beam modeling and hadron pro- 
duction at the target, uncertainties on neutrino cross 
sections, most of which are determined by in-situ cross- 
section measurements at MiniBooNE or other experimen- 
tal or theoretical sources, and uncertainties on detector 
modeling and reconstruction. By considering the varia- 
tion from each source of systematic uncertainty on the 9 e 
CCQE signal, background, and 9^ CCQE prediction as 
a function of E® E , a covariance matrix in bins of E® E 
is constructed, which includes correlations between 9 e 
CCQE (signal and background) and £7, CCQE. This co- 
variance matrix is used in the x 2 calculation of the oscil- 
lation fit. 

Figure [1] (top) shows the E® E distribution for 9 e 
CCQE observed data and background. A total of 144 
events pass the 9 e event selection requirements with 
200 < E® E < 3000 MeV. The data agree with the 
background prediction within systematic and statistical 
uncertainties. Fig. Q] (bottom) shows the event excess 
as a function of E~f . Also shown are expectations 
from the best £y, — > 9 e oscillation parameters returned 
by the fit and from two other sets of neutrino oscilla- 
tion parameters from the LSND allowed region The 
best oscillation fit for 200 < E$ E < 3000 MeV cor- 
responds to {Am 2 , sin 2 26>) = (4.42 cV 2 , 0.004), and 
has a x 2 of 18.2 for 16 degrees of freedom (DF), corre- 
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FIG. 2: The Q 2 (top panel) and cos(#) (bottom panel) dis- 
tributions for data (points with statistical errors) and back- 
grounds (histogram with constrained systematic errors) for 
E$ E > 200 MeV. Also shown are the expected distributions 
from intrinsic v e and v £ , and NC n° and A — > V7 back- 
grounds. 

sponding to a x 2 -probability of 31%. The null fit yields 
X 2 /DF = 24.5/18, with a x 2 -probability of 14%. A fit to 
475 < E$ E < 3000 MeV returns similar best-fit oscilla- 
tion parameters, (Am 2 , sin 2 29) = (4.42 eV 2 , 0.005), with 
X 2 /DF = 15.9/13 and a ^-probability of 25%. The null 
fit to 475 < E® E < 3000 MeV yields x 2 /DF = 22.2/15, 
with a x 2 -probability of 10%. The number of data, back- 
ground, and excess events for different E® E ranges are 
summarized in Table|TTJ No significant event excess is ob- 
served for Ef? E > 475 MeV. Furthermore, no significant 
excess is observed for E~} E < 475 MeV, to be compared 
to a 3. Oct excess observed for 200 < E$ E < 475 MeV in 
the v e appearance analysis 



TABLE II: The number of data, background, and excess events 
in the u e analysis for different E® E ranges. The correspond- 
ing numbers from the v £ analysis 0/ are on the right. The 
uncertainties include both statistical and constrained system- 
atic errors. 
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FIG. 3: Top: MiniBooNE 90% CL. limit (solid black) and 
sensitivity (dashed black) for events with E® E > 200 MeV, 
within a two neutrino z? M — ► v £ oscillation model. Also shown 
is the sensitivity for E® E > 475 MeV (dashed grey), and 
the limit from the BDT analysis (solid blue) for E$ E > 500 
MeV. Bottom: Limits from MiniBooNE for E® E > 200 MeV 
and E$ E > 475 MeV, KARMEN 0], and Bugey 0. The 
island contour in the bottom left corner is a lower sin 2 29 
limit from a fit to E® E > 475 MeV, excluding the points left 
of the line at 90% C.L.. The MiniBooNE and Bugey curves 
are 1-sided limits for sin 2 29 corresponding to A\ 2 = 1.64, 
while the KARMEN curve is a "unified approach" 2D contour. 
The shaded areas show the 90% and 99% CL. LSND allowed 
regions. 

The u e data also exhibit reasonable agreement with 
predicted background in other reconstructed kinematic 
variables. Fig. O shows the observed and predicted event 
distributions as functions of reconstructed Q 2 and cos(0) 
for 200 < E® E < 3000 MeV. Q 2 is determined from the 
energy of the outgoing lepton and its scattering angle 
with respect to the incident neutrino direction (9) assum- 
ing CCQE scattering. Also shown in the figures are the 
predicted distributions from NC ir° and A — ► V7 back- 
grounds, which arc events with a photon in the final state. 
The null x 2 values from these comparisons are both ac- 
ceptable, at x 2 /DF = 10.6/11 and x 2 /DF = 8.4/11 for 
Q 2 and cos(0), respectively. 

The absence of a significant excess allows MiniBooNE 



to place a limit on 



oscillations as shown in Fig. [3] 



The bottom panel of the figure shows the MiniBooNE 
limits obtained from fits to events with E® E > 200 MeV 
and E$ E > 475 MeV. Each 90% CL. limit on sin 2 29 is 
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obtained by a single-sided raster scan of the parameter 



Xb est fit 



< 1.64 cut is 



space, where a A% 2 

applied for each slice in Am 2 . The two limits are in 
agreement, with the one obtained for E® E > 200 MeV 
placing a stronger bound for low Am 2 oscillations, due to 
its slightly better sensitivity in that region (see top panel 
of Fig. [3]). At higher Am 2 values, both limits approach 
the corresponding sensitivities of the experiment, but at 
lower Am 2 both limits are noticeably worse due to the 
observed data fluctuation between 475 < E® E < 675 
MeV. The significance of that fluctuation in the 475 < 
E$ E < 675 MeV range is 2.8a (statistical © constrained 
systematic). 

Following 0], a secondary analysis based on Boosted 
Decision Trees (BDT) has been performed and used as a 
cross-check for the oscillation analysis in the energy re- 
gion E® E > 500 MeV, where the BDT analysis is not 
dominated by systematic uncertainties. No significant 
excess of events is observed with the BDT analysis, yield- 
ing the limit shown in the top panel of Fig. [3] Although 
the limit from the BDT analysis is not as stringent as 
the main result discussed above, the two analyses are 
complementary and yield consistent results. 

In summary, MiniBooNE observes no significant ex- 
cess of v e events in the energy region E® E > 200 MeV, 
for a data sample corresponding to 3.39 x 10 20 POT. 
Thus, with current statistics, MiniBooNE places a limit 
on two-neutrino — > v e oscillations shown by the black 
line in Fig. [3] The result is inconclusive with respect 
to small amplitude mixing suggested by the LSND data, 
but more antincutrino data, which are currently being 
collected, will provide additional information. More con- 
straints may also be provided b y th e off-axis NuMI beam 
data collected in MiniBooNE |22j. Interestingly, Mini- 
BooNE observes no significant excess of v e events in the 
low energy region 200 < E® E < 475 MeV. The absence 
of an excess at low energy in antincutrino mode should 
help distinguish between several hypotheses suggested as 
explanations for the low energy excess observed in neu- 
trino mode. 
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